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Abstract: In this paper, in the context of the quasi-magnetostatic approximation, we ex- 
amine incremental motions superimposed on a static finite deformation of a magnetoelastic 
material in the presence of an applied magnetic field. Explicit expressions are obtained 
for the associated magnetoacoustic (or magnetoelastic moduli) tensors in the case of an 
incompressible isotropic magnetoelastic material, and these are then used to study the 
propagation of incremental plane waves. The propagation condition is derived in terms of 
a generalized acoustic tensor and the results are illustrated by obtaining explicit formulas in 
two special cases: first, when the material is undeformed but subject to a uniform bias field 
and second for a prototype model of magnetoelastic interactions in the finite deformation 
regime. The results provide a basis for the experimental determination of the material pa- 
rameters of a magneto-sensitive elastomer from measurements of the speed of incremental 
waves for different pre-strains, bias magnetic fields, and directions of propagation. 
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1 INTRODUCTION 



This paper is in part motivated by the contributions of Michael Carroll to the subject of 
wave propagation effects in electroelastic and/or magnetoelastic materials, areas to which 
he made a range of contributions, including the papers [H 121 [3]. The subject has received 
new impetus in recent years because of the development of so-called 'smart materials', for 
which there is a strong interaction between mechanical and electric or magnetic responses, 
and, in particular, large deformations can be induced by an electric or a magnetic field, as, 
for example, in magneto-sensitive elastomers. The propagation of waves in such materials is 
therefore clearly of interest, but there is as yet relatively little literature that is concerned 
with electroacoustic or magnetoacoustic waves in finitely deformed and/or pre-stressed 
materials. Notable early exceptions are the paper by Yu and Tang which examined 
some limited aspects of waves in an infinite conducting magnetoelastic material subject to 
specific initial stresses, and two papers by De and Sengupta [5l |6] , in which surface and 
interfacial waves in an initially stressed conducting magnetoelastic material were studied. 
The formulation used in each case was rather specialized, however. 

A review of the state-of-the-art in the propagation of waves in magnetoelastic materials 
up to 1981 was provided by Maugin [7], in which paper the general equations governing 
magnetoelastic waves in deformable materials were developed, including, inter alia, the 
effects of finite deformation. Subsequent contributions have included the analysis of Abd- 
AUa and Maugin [8] , which used the general nonlinear equations of magnetoacoustics in the 
study of surface waves, the work of Boulanger ^ on inhomogeneous plane waves, and that 
of Hefni et al. (TU] on the nonlinear analysis of surface waves in an electrically conducting 
isotropic magnetoelastic half-space with a finite and constant magnetic field normal to the 
sagittal plane. 

With this background in mind a general formulation of the equations governing elec- 
tromechanical interactions in the finite deformation regime was reviewed by Ogden [TT] 
and applied briefly to the analysis of plane harmonic waves in a non-conducting magne- 
toelastic material within the quasi-magnetostatic approximation. This was in part based 
on and facilitated by the compact Lagrangian (or referential) formulation of the equations 
governing magnetoelastostatic interactions established by Dorfmann and Ogden [12]. In 
the static context the latter formulation was used in the construction of the constitutive 
laws and equations governing linearized incremental deformations superimposed on a finite 
deformation and an applied magnetic field (often referred to as a bias field) in a magne- 
toelastic material by Ottenio et al. [13] and employed in studying the combined effects 
of finite deformation and a magnetic field on the stability of a half-space. In the present 
paper we develop further the analysis in the paper by Ogden [TT] by providing a gen- 
eral three-dimensional analysis of incremental plane waves in an incompressible isotropic 
magnetoelastic material that is subject to both a finite deformation and a finite magnetic 
field. 

In Section |2]the basic equations governing static magnetoelastic interactions under finite 
deformations are summarized. This is followed in Section [3] by a derivation of the equa- 
tions governing incremental motions and the accompanying incremental magnetic field. 
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The incremental constitutive equations involve fourth-order, third-order and second-order 
magnetoelastic 'moduli' tensors that embody the material properties, and explicit expres- 
sions are given for these in the case of a general incompressible isotropic magnetoelastic 
material. 

Section H] is devoted to the study of incremental plane waves in an incompressible 
isotropic magnetoelastic material on the basis of the equations summarized in Section [31 
This is enabled by the introduction of a generalized acoustic tensor that governs the wave 
speed and the polarization of the mechanical displacement wave. Some algebraic manip- 
ulations enable us to find a symmetric representation of the acoustic tensor, leading to 
two possible transverse waves with mutually orthogonal polarizations. The possibility of 
combining these waves to form a circularly-polarized wave is also discussed. The mag- 
netic (induction) polarization vector is determined by a separate equation once the wave 
speed and mechanical polarization are found. These general results are then applied in 
two illustrative cases, for which the generalized acoustic tensor takes on relatively simple 
forms. We consider the propagation of plane waves, first in an undeformed material in 
the presence of a uniform bias magnetic field and, second in a homogeneously deformed 
material governed by a magnetoelastic constitutive law that is based on the Mooney-Rivlin 
material used in rubber elasticity, again in the presence of a uniform magnetic field. The 
latter constitutive law can be considered as a simple prototype model for describing the 
properties of magneto-sensitive elastomers. Explicit formulas for the wave speed and po- 
larization vectors are given and in some situations it is noted that the magnetic (induction) 
polarization vanishes, i.e. there is no disturbance to the underlying (bias) magnetic field 
by the mechanical wave. Section contains a brief concluding discussion. 

2 FINITE ELASTIC DEFORMATIONS IN THE 
PRESENCE OF A MAGNETOSTATIC FIELD 

We consider a magneto-elastic material, initially at rest in an unloaded and stress-free ref- 
erence configuration in the absence of an external magnetic field. We denote this configura- 
tion by Br- The material is subjected to the combined effects of a finite time-independent 
deformation, described by the deformation gradient tensor F, and a time-independent 
magnetic field vector H and associated magnetic induction vector B. The resulting cur- 
rent configuration is labelled B. We consider the material to be incompressible, so that 
det -F = 1, and we denote by p its mass density. In the equilibrium configuration B the 
total Cauchy stress tensor r satisfies the equation 

divr = (1) 

in the absence of mechanical body forces, while H and B satisfy the reduced Maxwell 
equations 

curl J/ = 0, divS = (2) 
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appropriate for magnetostatics, assuming that there are no free currents. In the above and 
in what follows div and curl denote the usual vector differential operators (with respect to 
B). 

By defining the total nominal stress hj T = F (for an incompressible material) and 
the Lagrangian counterparts of H and B by 

Hi = F^H, Bi = F 'B, (3) 

we may recast the equations ([T]) and ([2]) in Lagrangian form as 

DivT = 0, Cut\Hi = 0, DivBi = 0, (4) 

where Div and Curl are the Lagrangian counterparts of div and curl (with respect to Br), 
remembering that the incompressibility constraint det i*" = 1 is in force. 

There are many different ways in which the constitutive law for a magnetoelastic ma- 
terial may be formulated, but here we follow the formulation of Dorfmann and Ogden [12] , 
which is based on a total energy function, denoted by Q and defined per unit reference 
volume and expressed as a function of F and Bi: 

n = n{F,Bi). (5) 

Then, T and Hi are given by the simple formulas 

where p is a Lagrange multiplier necessitated by the incompressibility constraint. 

For an incompressible magnetoelastic material which is initially isotropic, it can be 
shown (see, for example, [12]) that f2 is a function of five independent invariants of the 
deformation and magnetic induction vector. The five invariants 

/i = trc, /2 = |[(trc)2-tr(c2)], h = Bi Bi, h = Bi cBi, I, = Bi-c'Bi (7) 

are typically used, where c = F'^F is the right Cauchy-Green deformation tensor. Then, 
T and Hi may be expanded in the forms 

T = -pF''^ + 2niF^ + 2n2{hF^-cF^) + '2^^5Bi®FBi + 2nG{Bi^FcBi + cBi®FBi), 

(8) 

and 

Hi = 2{n^Bi + n5cBi + nGc'Bi), (9) 

where Qk = dfl/dlk, k = 1,2, 4, 5, 6. Their Eulerian counterparts are 

r = -pi + 2Qib + 2fi2(/i6 - b^) + 2Q5B ®B + 2Qe{B ® bB + bB ® B), (10) 

and 

H = 2{yLib-^B + Vl^B + VLQbB), (11) 
where b = FF^ is the left Cauchy-Green tensor. 
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3 SUPERIMPOSED INCREMENTAL MOTIONS 

Superimposed on the static state described above we now consider an infinitesimal in- 
cremental motion for which the displacement, denoted u = u{x,t), is accompanied by 
an increment in the magnetic induction field and hence, via the constitutive equations 
(jH]) and, as appropriate, their isotropic specializations, induced increments in T and Hi. 
Here we are considering the quasi-magnetostatic approximation, so that time derivatives 
of the magnetic field and any electric field are neglected. Let a superposed dot indicate 
an increment. Then, the increments in F, T, Bi and Hi are denoted i^, T, Bi and Hi, 
respectively. We now refer these to the configuration B by the push-forward operations 
FP-^ = L = gradn, FT = Tq, FBi = Biq and F^Hi = Hiq, wherein the symbols 
with the attached subscript and the displacement gradient L are defined. 

The incremental forms u and Biq of the basic independent variables satisfy the in- 
cremental incompressibility condition and the incremental counterpart of ([2]) 2, which are, 
respectively, 

diva = 0, divBio = 0, (12) 

while, on taking the increments of (E]) and using the push-forward operations indicated 
above, we obtain the incremental constitutive laws 

fo = AoL + CoBio-pI + pL, Hio = C^L + CoBio, (13) 

where A^, Co and Co are, respectively, fourth-, third- and second-order tensors that depend 
on the material properties. Note that ^Aq is a generalization of the elasticity tensor in 
conventional elasticity theory, Co is a generalization of the inverse of the permeability 
tensor in magnetostatics, while Cq is a generalized magnetoelastic coupling tensor. The 
(Cartesian) components of these tensors are defined by 

(9^0 (92q /52Q 

where = {F^^)ai and the usual summation convention for repeated indices is used. 
We note the symmetries 

■A^piqj "^Oqjpi) ^ip\q Cpj|q, Coij Cpji. (-^^) 

Particular care is needed in using 'Cpi\q since it could also be written 'Cq\pi by interchanging 
the partial derivatives in f|T^ 9. The vertical bar is used to separate the pair of indices 
associated with the deformation gradient and the single index associated with the magnetic 
induction. See [13] for a derivation of these expressions in the purely static context, albeit 
in different notation. 

The incremental equation of motion is 

Divf = pu^tt, (16) 
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where is the material time derivative, and the incremental form of equation is 
Curli?/ = 0. When expressed in push- forward form on use of f|T3|) and the incompressibil- 
ity condition ( !T2|) these give 

dw{AoL + CqBio) - gradp + L^gi&dp = pu^u, curl(Co I, + Cq^^o) = 0. (17) 

For an isotropic material the components (fT4|) can be expressed in terms of the deriva- 
tives of fl with respect to the invariants ([7]). Using the results of Ottenio et al. [13] we 
find in turn 

A-opigj = 2 |r2i(5jj6pg + \2bipbjq — bigbjp — bijbpg + hSijbpg — 5jj(6^)pq] + Vt^dijBpBg 
+ ^6 [Sij{bB)pBq + 5ij{hB)qBp + bpgBiBj + bjpBiBq + bigBjBp + bijBpBq]} 
+ 4 {nnhpbjg + fi22(/i6 - b%p{hb - b^)jg + fii2 [bipihb - b^)jg + bjg{hb - b^),p] 
+ ^]l5 (bipBjB, + bjgB.Bp) + ^]25 [BiBp{hb - 6'),, + BjBg{hb - b%p\ 
+ n,5BiBjBpBg + ^]66 [ibB)iBp + (bS)pS,] [(6B),fi, + (6B),5,] 
+ ^]l6 [kpBg{bB)j + hpB,{bB)g + b,gBp{bB)i + bjgB,{bB)p] 
+ ^26 [(/i6 - b')ip [(^'S),5, + (6B),5,] + (/lb - b^) [(6B),5p + {bB)pB,]] 
+ [BiBpBg{bB)j + BiBpBj{bB)g + B.B^BpibB), + BjBgBi{bB)p]} , 

(18) 



Coiplq = 2 {05((5p5i?i + SiqBp) + ['^ig(b-B)p + 5pq{bB)i + fepq-Bj + 6ig-Bp]} 

+ 4 {niAp{b-^B)q + (]24(/i6 - 6')ip(b"'B), + n^5BiBpib-^B)q 
+ ^46 [(b-B),5p + {bB)pBi] {b-^B)q + r]i5&^p5, + Vl2^{hb - b%pBq 

+ n^^BiBpBq + [{bB)iBpBq + {bB)pBiBq + {bB)qBiBp] + VL^^bip{bB)q 

+ n^^ihb - b\p{bB)q + [{bB).,Bp + (bS)p5,] (6B) J , (19) 

and 

Coi, = 2 [(^4(&"')^i + ^^s^i, + + 4 {f]44(&"'-B)i(6-^B), + ^45 [{b-^B),B, 

+ (b"'B),-5i] + ^46 [{b-^BUbB)j + (b-^S),-(6B)i] + ^^r^B.B, 
+ ^56 [(&S)i5, + {bB)jBi] + fi66(bS)i(^'S)i} , (20) 

where Vtjk = d'^Vt/ dljdh, j, k G {1, 2, 4, 5, 6}. 



4 PLANE WAVES 

We now consider the basic (underlying) configuration to be uniform, consisting of a homo- 
geneous deformation and a uniform magnetic field, so that p is also uniform, the term in 
gradp in the equation of motion (fT7|l i vanishes, and the tensors A^, Cq, Cq are constants. 
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Moreover, we focus on the propagation of bulk infinitesimal homogeneous plane waves 
whose direction of propagation is defined by the unit vector n. We denote the wave speed 
by V. Thus, we seek solutions of the incremental equations f lT7|) in the form 

u = mf{n ■ X — vt), Bio = Q9{i^ " ^ — vt), p = P{n ■ x — vt), (21) 

where m and q are constant unit vectors in the directions of linearized polarizations for 
the mechanical and magnetic induction parts of the wave, respectively, and /, g, and P 
are single- variable functions of the argument n ■ x — vt. 

4.1 Magneto- acoustic tensors 

On use of equations (12T|) the equations (fT2|) yield 

m ■ n = 0, g ■ n = 0, (22) 

showing that both the incremental displacement and magnetic induction are transverse to 
the direction of propagation. 

The incremental equations f ll7l) can now be put in the form 

Q{n)mf" + R{n)qg' - P'n = pv^mf", n x [R{nfmf" + Coqg'] = 0, (23) 

where a prime signifies differentiation with respect to the considered argument and the 
components of Q{n), the acoustic tensor, and of R{ti), the magneto- acoustic tensor, are 
given by 

[Q{n)]-- = AopigjUpUg, [R{n)]-- = CoipijUp. (24) 

Note that, by f|T5l) . Q(n) is a symmetric second-order tensor, but in general R{n) is not 
symmetric. 

Now, by taking the dot product of (1231) i with n we obtain an expression for P', which 

is 

P' =[n- Q{n)m]f" + [n ■ R{n)q]g'. (25) 
Also, by noting that it follows from (123|) 9 that 

R{nfmf" + Coqg' = an, (26) 

for some scalar a, and taking the dot product of this with n we obtain 

a=[n-R{nfm]f" + {n-Coq)g'. (27) 

Substituting the expressions for P' and a back into equations (l23i) i and (126|) . respec- 
tively, we find 

I{n)Q{n)mf" + I{n)R{n)qg' = pv^mf", I{n)R^{n)mf" + I{n)Coqg' = 0, (28) 



7 



where I{ti) = / — n (g) n is the symmetric projection tensor on to the plane normal to n 
and is therefore the two-dimensional identity tensor in that plane. It has the properties 

I{nf = I{n), /(n)n = 0, I{n)m = m, I{n)q = q. (29) 

Using these properties we may rewrite (l28i) as 

Q{n)mf" + R{n)qg' = pv'mf", R{nfmf" + Co{n)qg' = 0, (30) 

where 

Q(n) = I{n)Q{n)I{n), R{n) = I (n) R{n)I (n) , Co(n) = /(n)Co/(n) (31) 

are the projections of the (three-dimensional) tensors Q{ti), R{n), Cq, respectively, on to 
the two-dimensional vector space normal to n. 

Let us assume that the incremental permeability tensor is positive definite and hence 
that its inverse Cq is also positive definite. It follows that Co{n) is positive definite on the 
two-dimensional space normal to n. Equation fl5Ul) 9 then yields 

qg' = -Co{n)-'R{nfmf", (32) 

and hence, on substitution of this into (|5U]) i . we obtain the equation 

r(n)m = pv'^m, (33) 

where 

f (n) = Q{n) - R{n)Co{n)~^R{nf, (34) 

which is a generalized acoustic tensor, sometimes referred to as the Green-Christoffel or 
simply Christoffel tensor. 

We remark that the Christoffel tensor T{n) in the eigenvalue problem fl33|) is symmetric, 
which ensures that the eigenvalues pv'^ are real and that the corresponding eigenvectors m 
are mutually orthogonal in the plane normal to n. And since r(n) is a two-dimensional 
operator there are at most two such eigenvalues. The eigenvalues are guaranteed to be 
positive (and hence the wave speeds real) if the inequality 

m-f(n)m>0, (35) 

holds for all unit vectors n, m such that n-m = 0. This is a generalization of the so-called 
strong ellipticity condition arising in pure elasticity theory, which, for an incompressible 
material, takes the form m ■ Q{n)m > 0. Note that since the eigenvalues pv"^ correspond 
to specific eigenvectors m for a given n their positiveness does not in general ensure that 
fl5S]) holds for all unit vectors m satisfying n • m = 0. 

By way of application of the above analysis we consider two examples. First, the 
propagation of waves in an undeformed magnetoelastic material in the presence of a finite 
applied magnetic field (often referred to as a bias field). Second, when there are both a 
finite deformation and a bias field present, we specialize the general analysis above to a 
prototype model that is appropriate for the description of magneto-sensitive elastomers. 
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4.2 Results for an undeformed material subject to a bias field 

When the material is undeformed [F = I) but subject to a uniform magnetic field the 
components of the moduli tensors are greatly simplified. For the tensor of elastic moduli 
we obtain 

Aopiqj = Oii6ip6jq + a26ij6pq + a^SiqSjp + a\ {6ipBjBq + 6jqBiBp) + a^SijBpBg 

+ ae{6pqBiBj + djpBiBq + 6iqBjBp) + arBiBjBpBq, (36) 

where the coefficients ai, . . . , 0:7 are constants given by 

ai =4(^11+4^22 + 4^12 + ^^2), ^2 = 2(^1 + ^2), a3 = -2fi2, as = 2(1^5 + 31]6), 

04 =4(^15 + 2f]25 + 2r]i6+ 41^26), ^6 = 2^6, 07 = 4(^55+4^56+4^66), (37) 

all the derivatives of Q being evaluated for F = I. Similarly, for the tensor Cq, we find 

Coipig = ll{SpqBi + 5iqBp) + 'y2BiBpBq + 'JsSipBq, (38) 



where 



and, for Co, 
where 



71 = 2(^5 + 2^]6), 72 = 4(r]45 + 2^]46 + ^55 + 3^56 + 2^66), 

73 = 4(^14 + 2Q24 + ^^15 + 2^25 + ^16 + 2^26), (39) 

Coij = Ci6ij + C2BiBj, (40) 



01=2(^4 + ^5 + ^6), 02 = 4(^44 + 2^45 + 2^46 + ^55+2^56 + ^66)- (41) 



Now, on substitution of these expressions into the definitions fl3T|) . we find that all three 
tensors Q{n), R{n), and Co(n) have the same structure, and are given by 

Q(n) = qiI+q2B®B, R{n) = (B ■ n) {-fj + -f2B ® B) , Co(n) = Cil+C2B®B, (42) 

where B = IB and 

gi = 2(^1 + ^2) + 2(^5 + 3^6) {B ■ nf , q2 = 2^6 + 4(^55 + 4^56 + 4^66) {B ■ nf . (43) 
We can easily derive the inverse of Co(n) as 

Co(n)-^ = -/ . ""'^ . B ® B, (44) 

^1 Cl [ci + C2B ■ Bj 

and then the Christoffel tensor is obtained as 

f{n) = aI + bB(^B, (45) 
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where 

7l/n .2 , I2 fry ^2 (Cl72 - C27l)^ .2 ..^n 

a = qi [B -n) , b = q2 [B ■ n) — [B ■ n) . (46) 

^1 ^2 CiC2 (^Ci + C2B ■ Bj 

This decomposition yields the two eigenvalues 

pt,2 = a, pv"^ = a + b B ■ B. (47) 

When these eigenvalues are substituted into the equation (l33l) we obtain the associated 
eigenvectors m. For pv'^ = a we obtain either S-m = or 6 = 0, while for pv"^ = a+b B B 
there are three possibilities: 6 = or S = or, provided B and B ■ m 0, m is 
aligned with B. These options are embraced by the following discussion. 

The first special case to consider is that of a wave travelling in a direction orthogonal 
to the magnetic induction vector, so that B ■ n = 0. Then, by (142|) 9. R{n) = 0, and it 
follows from (I30p 9 that g' = 0. Thus, in this case there is no incremental magnetic field, 
and the two wave speeds reduce to 

pv^ = 2{Qi + n2), pv^ = 2{ni + Q2) + 2 Qq B ■ B. (48) 

If -B 7^ these are distinct provided Qq ^ 0. li Qq = then circularly-polarized waves 
are possible, as will be discussed in the next section. Note that if B = then pv"^ = 
2{Qi + Q2) = p is a. repeated root corresponding to the classical elastic shear wave, where 
/i is the shear modulus of the material. 

The second special case concerns wave propagation in the direction of the applied 
magnetic induction vector, so that B = and n ■ B ^ 0. The two eigenvalues of the 
Christoffel tensor are then the same, leading to the conclusion that a circularly-polarized 
wave can always propagate in the direction of B in an undeformed magnetoelastic material. 
Moreover, g' ^ and the incremental mechanical motion is coupled with an incremental 
magnetic field, except in the special case for which + 2Qq = 0. 

A third special case can arise if 6 = 0, as indicated above, which can obtain under 
specific constitutive assumptions, an example of which is discussed in the following section: 
circularly-polarized waves can propagate in any direction in such a magnetoelastic material 
in the undeformed configuration. 

In the general case, for which n-By^O, nxBy^O, b^O, there are two linearly- 
polarized transverse waves, which may be combined to form an elliptically-polarized (but 
not circularly-polarized) wave. One is polarized in the direction of nx B and travels with 
speed [a/p]^/2j the other is polarized in the direction B = B — [B ■ n)n and travels with 
speed [{a + bB ■ B) / p]^/"^ . The existence of these waves is guaranteed if the strong ellipticity 
condition ( ^3]) holds. For the present specialization ( 15^ reduces to 

[r(n)m] ■ m = a + b{B ■ m)^ > for all unit vectors m, n such that m ■ n = 0, (49) 

where a and b depend on n through the product B ■ n. By choosing m so that B ■ m = 
and m = B /\B\ in turn it follows, in particular, that 

a > 0, a + bB ■ B > 0. (50) 
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In each case, g' ^ and the polarization of the incremental magnetic induction vector is 
aligned with that of the incremental mechanical wave, i.e. q = ±m. In this particular 
example it is easy to show that the inequalities fjSOjl imply the strong ellipticity inequalities 



4.3 Results for a Mooney-Rivlin magnetoelastic material 

For the second illustration of the general results we focus on a prototype form of the total 
energy function Q, which we refer to as the Mooney-Rivlin magnetoelastic material. This 
is given by 

n = \C{h - 3) + \D{h - 3) + \Kh + \i^h, (51) 

where C, D, k and z/ are constants. Using f|T8|) . f lTOl) . f l20|) . flMl) and fl3Tl) . specialized in 
respect of (EI]), we obtain the expressions 

Q{n) = [C{n ■ bn) + u{B ■ n^] I + D6"\ R{n) = u{B ■ n)I, Co(n) = nb'^ + ul, 

(52) 

where the hat denotes the projection of a tensor on to the two-dimensional vector space 
normal to n. 

In order to proceed we need the inverse of Co(n). This is obtained by using the two- 
dimensional Cayley-Hamilton theorem in the form 

Co(n) - [trCo(n)]/ + [det Co(n)]Co(n)-i = O, (53) 

which yields 

r r \~i v + Kii ^ K .... 

Z/^ -|- UKti + K^l2 V -|- VKlx -|- 

where i\ = tr(6~^) and ^2 = det(6~^). On substituting the expressions (l52l) into (l34l) we 
obtain the appropriate specialization of the generalized acoustic tensor as 

r(n) = Cn ■ (bn)I + Db'' + ^ '-^(^121 + i^h'^)- (55) 

Next, we take the dot product of equation fl33l) with k = n x m, where (m, k, n) 
forms an orthonormal basis. This yields the condition 

mb^k = mb^k = 0, (56) 

which means that m and k are along the principal axes of the elliptical section of the 
ellipsoid x ■ b~^x = 1 cut by the plane n ■ x = [S]. This is exactly the same situation 
as for the purely elastic Mooney-Rivlin material: two transverse waves may propagate in 
any direction n in the deformed solid, one polarized along m, the other along k. 

Finally, on taking the dot product of equation (l33|) with m and using (155|) gives the 
squared wave speed, say f^, via 



pvl = C{n . bn) + , , ''''''' {B ■ nf + 



11 



D 



u^K{B-nf 

+ vnii + ^2^2 



m ■ b m), (57) 



which is easily shown to be an eigenvalue of fl33p in the considered specialization. This 
expression can be simplified by noting that (156 p implies 



ii = m ■ b + k ■ b ^k = m ■ b + k ■ b ^k, 

22 = (^m ■ b'^rnj (k ■ b'^k^ = (m ■ b-^m) {k ■ b^k) , (5^ 



and hence we obtain 



pvi = Cin ■ bn) + D{m ■ b 'm) + '^^("^ ' ^ '"^) . . (59) 

u + n{m ■ b m) 

Here we note that we recover the result for the purely elastic case given by Boulanger and 
Hayes [H] not only when there is no applied magnetic induction {B = 0), but also when 
the wave propagates in the plane normal to the magnetic induction {B ■ n = 0). We also 
recover, by setting D = 0, the expression established by Ogden [TlJ for a neo-Hookean 
magnetoelastic material. Equation f l59|) provides a theoretical basis for determining the 
material constants by measuring the wave speed under different pre-strains and bias fields 
and hence an assessment of the magnetoacoustoelastic effects. 

We remark that Vk, the speed of the other wave propagating in the same direction, has 
a similar expression to f l59|) . with m replaced by k. Using (157|) . we find that the difference 
in the squared wave speeds is given by 



D H ^ 

+ vnii + K^i2 



{m b^m - k b 'k). (60) 



This expression provides an indication of the directions of the acoustic axes, along which 
circularly-polarized waves propagate. They correspond to a double root of the Christoffel 
tensor, when Vm = Vk, which is clearly equivalent to m ■ b^m = k ■ b^k. This condition 
means that the acoustic axes are along the normals to the planes of central circular section 
of the X ■ b^^x = 1 ellipsoid. Note that, for a general ellipsoid, amongst all the planes 
through its centre only two cut the ellipsoid in a circle; these are the planes of central 
circular section. The directions of the acoustic axes are thus dictated entirely by geometri- 
cal nonlinearities since they depend only on the pre-strain and not on any of the material 
parameters C, D, k or u. We refer to Boulanger and Hayes [13] for a detailed discussion 
of acoustic axes in the purely elastic case, which coincide with the acoustic axes in the 
present situation. See Figured] for an illustration. 

Note that in the undeformed specialization, b = I, and then 

pvl = pvl = C + D + {B ■ nf , (61) 

in agreement with the results of the previous section (since 6 = here). 

To complete the picture, we must determine the incremental magnetic induction, using 
(|32|). We obtain 



+ unii + K^i2)g'q = —iy{B ■ n) + Kii)m — nb f". (62) 
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Figure 1: Depiction of the circles of central circular section on the ellipsoid x ■ x = 1 
with the left Cauchy-Green deformation tensor b diagonal with entries in the proportion 
4:9:1. The acoustic axes, along which circularly-polarized bulk waves may propagate, 
are the normals to the two circles in the case of a Mooney-Rivlin magnetoelastic material, 
independently of the values of its four material parameters. 

In the special case for which the wave propagates in a direction normal to the magnetic 
induction, B ■ n = and hence g' = 0, i.e. there is no incremental magnetic induction 
field. Otherwise, we see from the above that q is parallel to m since, according to fl56l) . 

the only non-zero component of b~^m is along m. Thus, q = ±m, i.e. the mechanical 
and magnetic polarizations are aligned. 

5 CONCLUDING REMARKS 

In Section H] we have derived an explicit expression for the Green- Christoffel tensor gov- 
erning the propagation of small-amplitude magnetoelastic waves in a magneto-sensitive 
hyperelastic solid subject to a finite pre-strain and a magnetic bias field. Although the 
acoustic tensor for an incompressible material is not symmetric a priori, it is always possi- 
ble to obtain an equivalent, symmetric representation, by projecting all tensors on to the 
vector space normal to the direction of propagation. This has been established in general 
for internally constrained elastic solids by Scott and Hayes [12], and we were able to ex- 
tend the symmetrization process to the present context of magnetoelastic coupling. As a 
consequence, the resulting Christoffel tensor is two-dimensional, and the propagation of 
waves is governed by a symmetric second-order eigenvalue problem, which can always be 
solved explicitly and entirely. 

The two examples of the application of this eigenvalue problem presented here provide a 
theoretical basis for determining the material constants by measuring the wave speed under 
different pre-strains and bias fields and hence an assessment of the magnetoacoustoelastic 
effects. The following protocol can be employed. First, apply a magnetic bias but maintain 
the solid in an undeformed state; then, equation fH7|) gives a direct way of evaluating the 



13 



constants a and h. A revealing specialization occurs when the wave propagates along a 
normal to the applied magnetic induction vector. Then, the speeds of the two transverse 
waves should be distinct according to fH8|) unless VLq = 0. If they are equal, then the 
prototype fl5T]) is a candidate for the modelling of the material in question, and its material 
parameters C, D, n, v can be determined by using formula fl59|l under varying values for 
n, 6, and B. In particular, the directions of the acoustic axes should not change when 
\'B\ varies at fixed n, b, and Bj\'B\. Otherwise, a different, more sophisticated candidate 
must be considered in order to model the material properties. 
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